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Aims. Spectroscopic monitoring of BL Lac objects is a difficult task that nonetheless can provide important information on the 
different components of the active galactic nucleus. 

Methods. We performed optical spectroscopic monitoring of the BL Lac object AO 0235+164 [z = 0.94) with the VLT and TNG 
£^ \ telescopes from Aug. 2003 to Dec. 2004, during an extended WEBT campaign. The flux of this source is both contaminated 
I . and absorbed by a foreground galactic system at z = 0.524, the stars of which can act as gravitational micro-lenses. 

Results. In this period the object was in an optically faint, though variable state, and a broad Mgn emission line was visible at 
[* — . all epochs. The spectroscopic analysis reveals an overall variation in the Mgn line flux of a factor 1.9, while the corresponding 
f""*) ■ continuum flux density changed by a factor 4.3. Most likely, the photoionising radiation can be identified with the emission 
component that was earlier recognised to be present as a UV-soft-X-ray bump in the source spectral energy distribution and 
that is visible in the optical domain only in very faint optical states. We estimate an upper limit to the broad line region (BLR) 
size of a few light months from the historical minimum brightness level; from this we infer the maximum amplification of the 
O ' Mgn line predicted by the microlensing scenario. 

Conclusions. Unless we have strongly overestimated the size of the BLR, only very massive stars could significantly magnify 
^2 the broad Mgn emission line, but the time scale of variations due to these (rare) events would be of several years. In contrast, 
. . \ the continuum flux, coming from much smaller emission regions in the jet, could be affected by microlensing from the more 
■ plausible MACHO deflectors, with variability time scales of the order of some months. 
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' - galaxies: jets - galaxies: quasars: general 
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1. Introduction 

The active galactic nuclei (AGNs) known as BL Lac ob- 
jects show extreme variability at all wavelengths from 
the radio band to the 7-rays. The commonly accepted 
paradigm foresees that their non-thermal emission comes 
from a relativistic plasma jet pointin g towards the ob- 
serve r with a small viewing ang le (|Blandford k Reesl . 
19781) . The consequent relativistic Doppler effect would 
be responsible for flux enhancement and contraction of 
the intrinsic variability time scales. 



to: 



C. M. Raiteri, e-mail: 



Send offprint requests 
raiteri@to . astro . it 

* Based on observations collected at the European Southern 
Observatory, Chile (ESO Programme 71.A-0174), and on ob- 
servations made with the Italian Telescopio Nazionale Galileo 
(TNG) operated on the island of La Palma by the Fundacion 
Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) 
at the Spanish Observatorio del Roque de los Muchachos of 
the Instituto de Astrofisica de Canarias 



These objects by definition show no emission lines or 
only weak ones, w ith equivalent width s not exceeding 5 A 



in the rest frame (jStickel et al 
emission lines 



19911 ) . However, stronger 
in particular broad ones, have occasion- 



ally been dete cted in the spectra of a few BL Lacs (e.g. 
AO 0235+164 ICohen et al.lll987t iNilsson et alJll99rj BL 
Lacer tae itself, IVermeulen et al.lll995l ; ICorbett et al.lll996t 
2000) in optically faint states, thereby putting their clas- 



sification into question. 

Those broad emission lines recognisable in the spectra 
of AGN are produced in the so-called "broad line region" 
(BLR), which is believed to be photoionised by thermal 
radiation from the accretion disc. In the case of BL Lac 
objects, the possibility that photoionisation can be due to 
the beamed radiatio n from the jet is still under debate 
(ICorbett et all l2000h . 

Because of the dramatic flux vari ability and absence of 
strong lines. lOstriker fc Vietri (1985) suggested microlens- 
ing as a possible explanation of the BL Lac phenomenon. 
According to this interpretation, the moving stars of a 
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galaxy on the line of sight would act as gravitational 
microlenses, producing amplification and the consequent 
variation in the flux emitted from a compact underlying 
source. A microlensing-induced flare would be observed 
simultaneously at different wavelengths with time scales 
that depend on the size of the regions emitting at the 
various frequencies. On the contrary, no amplification of 
the lines is expected, if the size of the BLR is of the or- 
der of 0.1-1 pc, as classically assumed. Actually, it was 
found that the BLR size scales with the continuum lu- 
minosity and that for some objects it can be rather com 



pact, of the order of som e light days-weeks (jWandel et al- 
199alKaspi et al.l . l2000h . In this case, microlensing would 



also produce noticeabl e variations in the emission line 
strengths and profiles (jAbaias et all [2002). One would 



thus observe these changes at the same time as the con- 
tinuum ones and would find that the line flux correlates 
with the continuum flux. It is worth mentioning that the 
observed difference between continuum and broad emis- 
sion line flux ratios of lensed quasar images has been in- 
terpreted as due to microlensing, le ading to an estimate 



of the BLR s i ze down to 9 light davs flMetcalf et all [2004; 



Wavth et all l2005t IKeeton et all . [2006) 



In this context, the BL Lac object AO 0235+164 at 
z = 0.94 is very interesting. It is well known for its strong 
flux variability at all wavelengths. In the optical band 
it has shown brightness variations up to ~ 5 mag in 6 
months, while radio flux changes up to a factor ~ 75 
in less than 1 year have been observed. In some cases 
simultaneous opt i cal and radio outbursts were observed 
dBalonek k Pent! . Il98d IWebb et all l200d IRaiteri et af 



2001), with a possibl e quasi-period i city o f the major radio 



(and optical) events (jRaiteri et all 12001V 



Some observational evidence can be reconciled with 
the intrinsic variability scenario for this source only if 
very extreme relativistic conditions are assumed: Doppler 



pers ( 


7 ujisawa ct al., 


1999 


2000; 


Jorstad et all 


200111 



Kraus et al. , 1999; iFrey et al 



Interstellar scattering would 
be a viable extrinsic explanation for the intraday radio 
variability, at least in some cases, but it cannot explain 
variability at the higher frequencies. The alternative ex- 
trinsic scenario, invoking microlensing, appears attrac- 
tive, since optical imaging and spectroscopic observations 
have revealed foregrou nd absorbing s ystems at z = 0.52 4 
and z = 0.851 (see e.g. Cohen et all Il987l ; iNilsson et al 



19961 ). Absorption in excess of Galactic absorption is in 
deed re quired to fit th e X-ray spectra of this source 



(sec eg 



Madeiski et al.. 19961 : Junkkarinen et al. . 2004 : 



Raiteri et all l2006ri ). The environment of AO 0235+164 



is rather complex (see Fig. [1}: the z — 0.524 absorber is 
probably an elliptical galaxy lying l'/3 to the eas10, which 
is very likely interacting with a Seyfert galaxy 2" south 
of the source (G2). This southern AGN, named ELISA 



1 Following the notation bv lNilsson et a.1.1 1 19961 . it was indi- 
cated as Gl in the figure, but it is not discernible from the 



Fig. 1. Detail of the i?-band frame of the AO 0235+164 
field taken with a 15 s exposure at VLT. The image is l'x 
V centred on the source; north is up and ea s t is on the 
left. Following the notation by Nilsson et al. ( 19961 ). Gl 
to G4 are the four galaxies at z = 0.524; the intervening 
galaxy, Gl, is not resolved from the source; the southern 
AGN, ELISA (G2), is instead clearly visible. 



by IRaiteri et al. ( 20051 ) , can noticeably affect the source 
photometry during faint states, especially in the blue part 
of the optical spectrum. The f our galaxies G1-G4 belong 
to the same z = 0.524 group (|Stickel et allll993h . which 
might include other close objects with still unknown red- 
shifts. 

The possibility that microlensing affects the AO 
0235+164 e mission was the subject of a number of works 
in the past dKayserlll988tlStickel et al 1 1 19881: [Saustl. 1 1 992 



Abraham et all Il993t IKraus et all 1 19991 IWebb et al. . 
20001) . which each faced the problem with different tech- 
niques (mainly imaging and multiwavelength photometric 
monitoring). The results of these stu dies are contradic- 
tory. According to Stickel et al. ( 19881 ). the presence of a 
foreground e astern galaxy very close to the line of sight, 
the fact that ICohen et alj ( 1987 ) observed that only mild 
line flux variations correspond to strong continuum-flux 



densit y changes, and the observation by iBalonek k Dent 
(|l980h of contemporaneous optical and radio outbursts, 
all support a microlensing origin of the violent variability 
of this source. Also the new VLBI components emitted af- 
ter outbursts are interpret ed as images o f the lensed core 
of the source. In contrast, Kavserl (|l988l ) showed how the 



microlensing scenario is very unlikely on the basis of sev- 
eral observational data. In particular, he calculated that 
the separation of the new VLBI components from the core 

f . 



would imply a non-stellar deflector m ass of 5.7 x 10 5 M, 



According to lAbraham et al. ( 1993f) . microlensing does 
not appea r to be a likely expla nation for the source proper- 
ties, while Kraus et al. ( 19991 ) claimed that it can account 
for the peculiar variations they observed at three radio 
frequencies and in the opti cal band du r ing a three-week 
campaign in 1992. Finally, Webb et al.l (2000) suggested 
that the events in 1975 and 1997, which were characterised 
by simultaneous optical and radio outbursts, might be as- 
cribed to microlensing, while others are more likely due to 
intrinsic causes. 

Actually, a careful i nspection of the hist orical optical 
and radio light curves (jRaiteri et all l2006af ) reveals that 
simultaneity, as well as symmetry of events, are features 
difficult to verify, either because of poor sampling (espe- 
cially lack of optical data on long time intervals due to 
solar conjunctions), or due to the possible overlap of dif- 
ferent events. Even for the best-sampled outburst in 1997- 
1998, we can only say that there was a strong flux increase 
in both the optical and radio bands, but the details of the 
outburst (exact time of the peak and shape at the various 
frequencies) largely remain unknown. 
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Investigation of the microlensing scenario through 
a comparison between line and con tinuum flux varia- 
tions was performed by ISaust (1992), who observed AO 
0235+164 in both the V band and a narrow-band filter 
centred on the Mgn line. He found that the emission- 
line variability was comparable to the continuum one and 
that it was stronger than expected from microlen sing. This 
conclu sion, however, was based on the models of lNemirofl 
(1 19881). which assu med large BLR sizes. More recently, 



Abaias et al.l (|2002i ) has shown that the global amplifica- 



tion of broad emission lines induced by microlensing events 
could be noticeabl e, if the BLR has a size o f a few light 
days, as found by IWandel et all (|l999h and iKaspi et al 
( 2000) for several AGNs. 



In this paper we present an analysis of AO 0235+164 
optical spectra taken with the VLT and TNG telescopes 
in 2003-2004, when the source was in a faint state. In this 
period, a huge multiwavelength observing effort was car- 
ried out by the Whole Earth Blazar Telescope (WEBTjl 
collaboration to continuously monitor the source emission 
variability in the optical and radio bands. Furthermore, 
three pointings of the XMM-Newton satellite provided 
information in the X-ray and UV bands, allowing us 
to construct the broad-band spectr al energy distribu- 



tions (S EDs) with simultaneous data ( Raiteri et al. . 20051 



tions (SI 
2006bllal ) 



The analysis of multiepoch SEDs suggested the 
existence of an extra component besides the low-energy 
synchrotron one and the high-energy inverse-Compton 
one. This extra component, peaking in the far-UV-soft- 
X-ray domain, can be explained in terms of an additional 
synchrotron component coming from an inner part of the 
jet with respect to that emitting the observed radio-to- 
optical radiation. However, the hypothesis of a thermal 
component coming from an accretion disc cannot be ruled 
out, particularly because of the possible detection of a 
strongly redshifted Fe Ka line in some of the X-ray spec- 
tra, which should b e produced in the in ner part of the 
disc, as discussed bv lRaiteri et al. ( 2006bl ). 



2. Observations and data reduction 

Optical spectra were taken in the period August 2003 - 
December 2004 at the 8 m diameter Very Large Telescope 
(VLT) and in service mode at the 3.6 m Telescopio 
Nazionale Galileo (TNG) . The logbook of the observations 
is shown in Table[T] where Column 1 reports the telescope, 
Column 2 the instrument, Column 3 the observing date, 
Column 4 the value of seeing derived from the acquisition 
frame(s), and Column 4 the total integration time. This 
was obtained in one exposure at VLT, while three sub- 
sequent exposures were acquired at TNG, in order to be 
able to get rid of cosmic rays by combining the images. 

All spectra were obtained with a 1" slit aligned along 
the north-south direction, so that the light of both the 
source and ELISA entered the spectrograph. 



The data reduction was carried out with the IRAI0 
package, following a standard procedure. Spectra were 
bias-subtracted and corrected with dome flat-field frames. 
Arc-lamp exposures were used for wavelength calibration; 
two-dimensional mapping of pixels into wavelengths was 
derived by running the routines identify, reidentify, 
and fitcoords, which allows correction for geomet- 
rical distortions. Wavelength-calibrated source spectra 
were then obtained with the transform routine. Sky- 
subtracted, one-dimensional source spectra were extracted 
with apall. A small aperture was used in order to avoid 
the inclusion of ELISA. From the VLT frame we could 
also extract the ELISA spectrum. 

Spectra of spectrophotometric standard stars were 
used for flux calibration, which was performed with the 
routines standard, sensfunc, and calibrate. 

We checked the correctness of flux calibration using 
field photometry. Indeed, for all observations, acquisition 
frames in at least the R band were taken, so that differen- 
tial photometry was performed with the routine phot with 
respect to comparison stars in the same field of the source. 
The high space resolution (0'.'25 per pixel at VLT and 
0'.'275 per pixel at TNG), together with the relatively small 
seeing values (see Table [T]), allowed us to avoid the inclu- 
sion of the ELISA contribution by using an aperture ra- 
dius of about 1". This is an important point, since ELISA 
can significantly affect the photometry of AO 0235+164, 
especially when th e source is in a faint state (see e.g. 
Raiteri et al. I, 120051) . For the magnitude calibration, we 
used the photometric sequences bv lSmith et al.l (|l985l ) and 
Fiorucci et al.l (|1998l ); the values of the source magnitude 
obtained with respect to the brighest, non-saturated com- 
parison stars were averaged, and the uncertainty reported 
in Table 4 takes into account both the standard deviation 
of these values from the mean and the instrumental error. 

The results of the photometric measures are reported 
in Tabled] As one can see, the source was always found in 
faint states, but at different levels, ranging from R ~ 17.5 
on November 15, 2003, to R ~ 18.8 on December 6, 2004. 
Figure [2] shows the i?-band light curve of AO 0235+164 
obtai ned during the 2003- 2005 WEBT campaign (adapted 
from iRaiteri et ai1l2006al) . 

To check the correctness of the spectroscopic flux cali- 
bration, we transformed the i?-band magnitudes into flux 
densities and compared them with the spectroscopic con- 
tinua around 6410 A. In order to verify that this proce- 
dure also leads to a good calibration around the Mgn 
line, we derived V-band magnitudes either from V-band 
frames, when available, or though the mean c olour index 
V — R = 0.76 estimated by Raiteri et al.l (|2005l ). obtaining 
flux densities at ~ 5450 A. In a couple of cases, frames in 
other bands were also available, allowing for a more de- 



2 http://www. to.astro.it/blazars/webt/ (see e.g. 
Ivillata et al.l , l2004bl ,lal. 120061 ) 



3 IRAF is distributed by the National Optical Astronomy 
Observatories operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation 
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Table 1. Observation log. 
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Telescope 


Instrument 


Date 


Seeing ["] 


Exp. [ 


VLT4 


FORS2, grism 300V 


Aug. 28, 2003 


0.8 


900 


TNG 


DOLORES, grism LR-B 


Nov. 15, 2003 


1.3 


1800 


TNG 


DOLORES, grism LR-B 


Jan. 17, 2004 


1.3 


1800 


TNG 


DOLORES, grism LR-B 


Oct. 5, 2004 


1.2 


3600 


TNG 


DOLORES, grism LR-B 


Dec. 6, 2004 


1.4 


3600 



2003.5 



2004.0 



2004.5 



2005.0 



16 



17 



a; 18 



f 11 In 




19 



20 r 

2800 2900 




3000 3100 3200 

Julian Date - 2450000 



3300 



3400 



Fig. 2. The 7?-band light curve of AO 
0235+164 during the 2003 -2005 WEBT 
campa ign (adapted from iRaiteri et al.l 
l2006af ). The dates of the spectroscopic 
observations, as well as the correspond- 
ing ii-band values derived from the ac- 
quisition images, have been indicated as 
vertical dotted lines and blue crosses. 



tailed comparison between spectroscopic and photometric 
results. 



2.1. VLT spectra 

The VLT data were acquired on August 28, 2003 under 
photometric conditions with sub-arcsecond seeing. Grism 
300V was used, with 11 A dispersion and 5900 A central 
wavelength. The CCD scale was 0'.'25 per pixel, so that the 
spectrum of ELISA, the southern AGN located 2" south, 
is well-resolved from the source one. Flux calibration was 
performed with the spectroscopic standard star LTT 1020, 
which was also used to remove the prominent telluric ab- 
sorption feature at 7550-7700 A. 

Figure [3] shows the calibrated source spectrum, where 
the continuum is traced well from about 3800 A to about 
9200 A. The spectrum shows a number of clear emission 
features at both z = 0.94 and z = 0.524 (see Table |2J). In 
particular, we can distinguish a broad Mgn emission line 
at A ~ 5430 A, together with other narrow lines emitted 
from the source: [Nev], [Oil], Nem, while the identifica- 
tion of the H<5 and H7 lines appears more uncertain. The 
full width half maximum (FWHM) of the Mg 11 emission 
line is ~ 3500 k m s" 1 , a bit la r ger th an the value of 3100 
kms -1 found by ICohen et al.l ( 1987t ); on the other hand, 



Table 2. Emission lines in the AO 0235+164 spectrum. 



Line 


Aobs 


z 


Flux 


FWHM 




[A] 


[10- 


16 erg cm -2 s _1 ] 


[kms" 1 ] 


A2800 Mgn 


5431 


0.940 


14 


3500 


A3426 [Nev] 


6640 


0.938 


1.2 


700 


A3727 [0 11] 


7225 


0.939 


2 


600 


A3869 [Nem] 


7503 


0.939 


2 


700 


A4102 m 


7958 


0.940 


4 


3200 


A4340 H7 


8439 


0.944 


9 


3200 


A3727 [O n] 


5682 


0.524 


1.3 


800 


A4861 H/3 


7407 


0.524 


1.8 


900 


A5007 [O m] 


7631 


0.524 


1.8 


900 



the 3200 kms 1 velocities we found for the US and H7 
lines are sma l ler bu t consistent with those measured by 
Nilsson et ail (|l996l ). 

The emission lines of [0 11] and H/3 at z = 0.524 belong 



to the foreground galaxy located at about l'/3 east of the 
source. Two pronounced features, probably due to Fe 11 
and Mgn absorption from the same galaxy at z = 0.524 
are evident on the red side of the spectrum, while on the 
blue side one can recognise Na 1 D absorption at z = 0.524. 
Moreover, t he absorption featu re at 5183 A was already 
ascribed by I Cohen etail (|l987l ) to Mgn at z = 0.851. 
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Fig. 3. The VLT spectrum of AO 
0235+164 acquired on August 28, 2003. 
Several emission features at z = 0.94, 
as well as z — 0.524, are evident, plus 
absorption lines at z — 0.524 and 2 = 
0.851. Telluric absorption bands are in- 
dicated by a "T" label, and the 7550- 
7700 A telluric feature has been cor- 
rected for. 



Table 3. Main emission lines in the ELISA spectrum. 2.2. TNG spectra 



Line 


Aobs 


Flux 


FWHM 




[A] 


[10" 17 erg cm^s" 1 ] 


[kms" 1 ] 


A2800 Mgll 


4260 


80 


4000 


A3426 [Nev] 


5220 


12 


800 


A3727 [0 11] 


5682 


35 


800 


A3869 [Nem] 


5895 


5 


800 


A4102 m 


6252 


5 




A4340 H7 


6615 


8 




A4861 H/3 n 


7410 


15 


700 


A4861 H/3 b 


7410 


25 


3000 


A4959 [Om] 


7556 


25 


700 


A5007 [O in] 


7630 


80 


700 


A5199 [Ni] 


7915 


4 


600 


A6300 [Oi] 


9606 


7 


700 



The spectrum of ELISA is shown in Fig.[4j Both broad 
and narrow emission lines are clearly recognisable. In par- 
tic ular, one can s ee the H/3 and [O Hi] lines already n oticed 
bv ICohen etaD (jl987h and bv INilsson et all (|l996h . The 
latter authors also identified the H7 line. 

All the other features highlighted in Fig.0]are reported 
for the first time here; the strong broad Mgll and the nar- 
row [0 11] lines are the most striking features. Table[3]sum- 
marises the identified emission lines in the ELISA spec- 
trum, together with the derived flux. A double Gaussian fit 
to the H/3 line to measure both the narrow and the broad 
co mponents yield s resul ts in agreement with those derived 
bv INilsson et all (jl996l) . We notice that the FWHM o f 



[0 11] and A5007 [O ill] measured by ICohen et alJ (jl987t ) 
were 300 and 250 kms" 1 , r espectively, and the F WHM of 



the [O III] lines measured bv lNilsson et al. (1996) were 400 
kms" 1 , somewhat smaller than the ones derived from our 
spectrum. 



The TNG spectra were obtained with DOLORES and 
grism LR-B. The resolution is 11 A. The scale of 
DOLORES is 0'.'275 per pixel, so that the source and 
ELISA spectra are resolved also in the TNG frames. 
However, the ELISA spectrum cannot be extracted be- 
cause of the low S/N ratio. The total integration time 
(reported in Table [I} of each observation was obtained 
with three subsequent exposures. The three corresponding 
spectra were reduced and calibrated separately, and then 
summed up in order to increase the S/N ratio. The mea- 
surements of the Mgll equivalent width (EW) and flux 
were performed on each single exposure; the values re- 
ported in Table [4] are the average values obtained over the 
exposures, and the uncertainty corresponds to the stan- 
dard deviation. These mean values are in fair agreement 
with those obtained by measuring these quantities directly 
on the summed spectra. The final summed TNG spectra 
are shown in Fig. [5l 

Sky subtraction on October 5 and December 6, 2004 
was poor because of CCD charge-dragging problems. In 
these two epochs a complete sequence of UBVRI images 
was taken, so that it was possible to check the flux cali- 
bration at different wavelengths in order to be sure that 
the flux around the Mgll line was correct. 

Figure [5] also shows flux densities derived from 
UBVRI photometry for these last two epochs. One can 
notice the deviation from a linear trend in the bluer part 
of the spectrum. This deviation is also confirmed by the 
shape of the spectroscopic continuum, even though the 
photometric data on December 6 give a higher flux den- 
sity in U and B bands than do the spectroscopic data. 
In any case, the data suggest a blue excess on December 
6, when the source was in the faintest state. This excess 
would be amplified when correcting for both Galactic and 
foreground galaxy absorption, and is likely to be the sig- 
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Fig. 4. The VLT spectrum of ELISA 
acquired on August 28, 2003. Several 
emission features are evident; telluric 
absorption bands are indicated by a 
"T" label; the 7550-7700 A telluric fea- 
ture has been corrected for. 



Table 4. i?-band magnitude, Mgn equivalent width and flux, corresponding continuum flux density, and FWHM of 
the Mgn line derived from the spectra presented in this paper. EWs and FWHMs are in the observed frame. 



Observation 


R 
[mag] 


Mgii EW 

[A] 


Mg ii Flux 
[10~ 16 erg cm" 2 s _1 ] 


Cont. Flux Density 
[10" 16 erg cm" 2 s _1 A" 1 ] 


FWHM 
[kms^ 1 ] 


VLT, Aug. 28, 2003 


17.90 ± 0.05 


10 ±1 


14 ± 1 


1.4 ±0.1 


3500 


TNG, Nov. 15, 2003 


17.49 ±0.03 


7.1 ±0.9 


13.8 ± 1.6 


1.94 ±0.05 


3300 


TNG, Jan. 17, 2004 


18.22 ±0.06 


12.5 ± 1.3 


11.3 ± 1.1 


0.90 ± 0.06 


3300 


TNG, Oct. 5, 2004 


17.93 ±0.02 


9.5 ±0.8 


11.1 ±0.9 


1.17 ± 0.02 


3100 


TNG, Dec. 6, 2004 


18.79 ±0.02 


15.7 ± 1.2 


7.2 ±0.5 


0.45 ± 0.02 


3400 



nature of the low-energy tail of the e xtra component sug- 
gested bv lRaiteri et"alT(|2005l l2006bl lah , which has already 
been mentioned in the introduction. 



3. Ionising source 

Table [4] summarises the results of the spectroscopic anal- 
ysis on the Mgn line. For each spectrum it reports: the 
i?-band magnitude, the Mgn line equivalent width and 
flux, the corresponding continuum flux density, and the 
line FWHM. One can see that the continuum flux den- 
sity varied by a factor 4.3, while the Mgn flux changed 
by a factor 1.9. The results are also presented in Fig. [6j 
which shows the Mgn flux as a function of the contin- 
uum flux density, the Mg n equivalent width as a function 
of the continuum flux density, and the Mgn equivalent 
wi dth versus the J?-ba nd magnitude. The values obtained 
bv lCohen et~aT1 (|l987h are also show n for a compa r ison. 
The main result, confirming the ICohen et al.l (J1987) 



claim, is that the Mgn flux is variable. However, its rela- 
tionship with the continuum flux density, if any, is not sim- 
ple. The Mgii equivalent width seems to depend linearly 
on the i?-band magnitude; the grey lines in the various 
panels represent this linear relationship: EW = a ± bR, 
with best-fit values a = -112.76 and b = 6.84. The R- 



band magnitude has been converted into the continuum 
flux density around the Mgn line, i.e. in the V band, 
through the V — R index. The three lines in the first and 
second panels correspond to the ave rage colour index an d 
to 1 a deviations from it according to lRaiteri et al. (2005): 
V — R = 0.76 ±0.10. For a comparison, in the second 
panel we also plotted the case where the Mgn line flux 
is assumed to be constant (EW cx F^ t ) and equal to 
its average value (10.94 x 10 _16 ergcm -2 s -1 ; red dashed 
line). 

From Fig. [6] one may infer that the line flux increases 
with the continuum flux density until a certain brightness 
level, after which it cannot grow further and may even de- 
crease. We notice that this first increasing trend depends 
on just one point, but further observations would be re- 
quired to confirm it. 

Corbett et ahl (|2000h discussed the possible BLR pho- 



toionisation source in BL Lacertae. They analysed eight 
optical spectra taken in 1995-1997 and measured the Ha 
broad emission line. The behaviour of its equivalent width 
as a function of the continuum flux density was well-fitted 
by an F c ~* t law. The interpretation of this lack of response 
by the line strength to the continuum variations in terms 
of photoionisation by Doppler-boosted synchrotron con- 
tinuum requires a high covering fraction of about 35%. 
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Fig. 5. The TNG spectra of AO 
0235+164 acquired on November 15, 

2003 (black), January 17, 2004 (green), 
October 5, 2004 (red), and December 6, 

2004 (blue). Flux densities derived from 
UBVRI photometry in the last two 
epochs are shown to compare the spec- 
troscopic with the photometric shape of 
the spectra. 
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Fig. 6. Results of the AO 0235+164 
spectroscopic monitoring at the VLT 
and TNG: Mgn line flux or equivalent 
width plotted against continuum flux 
density o r fl-band magnitud e; the re- 
sults by I Cohen et alj (1987) are also 
plotted for comparison. See text for fur- 
ther details. 



Although this scenario cannot be ruled out, the authors 
favoured photoionisation by a second continuum source 
not observable in the optical band. In the case that it 
is a hot accretion disc, they estimated a temperature 
T > 1.2 x 10 5 K. 

For AO 0235+164 this second continuum source can be 
identified with the UV- s oft-X-r ay bump in the SED dis- 
cussed by iRaiteri et al.l (|2006al lbl) . From a visual inspec- 
tion of the SEDs corresponding to different epochs, one 
may infer that this component is variable, but its varia- 
tions do not appear to be necessarily correlated with the 
optical ones. 

Thus, for AO 0235+164 we can envisage two scenar- 
ios: i) the BLR is photoionised by the jet emission; the 
line flux grows with the continuum flux but only until a 



kind of saturation level is reached, perhaps corresponding 
to complete ionisation; ii) the photoionisation source is 
another component, producing the UV-soft-X-ray bump 
in the SED; it can be detected in the optical domain when 
the lower-energy synchrotron component is weak, i.e. dur- 
ing very faint optical states. In this case the line flux is 
independent of the continuum flux, since they have differ- 
ent origins, and that we see a weaker line in the faintest 
state is just by chance. This component can be either 
emission from a thermal disc or an additional synchrotron 
component coming from an inner part of the jet with re- 
spect to the one producing the radio-optical emission. The 
thermal-disc option may be supported by the possible de- 
tection of a stro ngly redshifted Fe Kg line in the source 



X-ray spectrum ( Raiteri et al. . 2006bJ ) 
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Between these scenarios, we think that the second one, 
i.e. photoionisation by a different component with respect 
to what produces the optical continuum, is more likely, 
since there is observational evidence for its existence. 

4. Size of the BLR 

The size of the BLR in AGNs can be estimated 
with reverberation mapping techniques, but relation- 
ships have also been derived bet ween the BLR size , 

" 19991 : 



they used the equation derived by iMcLure fc Dunlopl 
(|20041 ) that applies the virial black hole mass estimate 



-Rblr, and the con t inuum luminosity (|Wandel et all „ 
Kaspi et all. l200rj IMcLure fc Dunlopl . |2004|). Inparticu- 



lar. lMcLure fc Dunlod (|2004l ) analysed the dependence of 
Rblr on the continuum luminosity at 3000 A in the source 
rest frame for high redshift objects (z > 0.7). They found 
that 

i?BLR = (18.5 ± 6.6)[A J L 300 o/10 37 W](°- 62±014 ) light days, 

where the continuum is supposed to be the ionising radi- 
ation of the BLR and can be estimated around the Mgn 
line. 

As discussed in the previous section, the most plausible 
photoionisation source of the BLR in AO 0235+164 is an 
extra component, which is visible in the optical frequency 
range when the source is faint. Indeed, we recognised a 
sort of blue excess in the faintest spectrum we analysed 
(taken on December 6, 2004). Hence, if we consider the 
lowest state of AO 0235+164 ever observed, it is plausible 
that a substantial fraction of the continuum there comes 
from the ionising radiation. 

The minimum brightness level of AO 0235+164 in the 

V band was observed in July 1996: V = 20.29 ± 0.170. 
Simultaneous data in other bands are not available, but 
there are i?-band data the day before and the day after, 
indicating a decrease in the source brightness. Hence, we 
can assume that their average value, R = 19.89 ± 0.10, is 
not far from the actual R value at the time the y-band 
datum was taken. In this case, one gets a colour index 

V — R — 0.40, which implies a rather flatter spectrum than 
usually found for this source. This can be an indication 
that the higher-energy component was dominant. 

If we consider absorption by both the Galaxy and 
the foregrou nd absorber at z = 0.524 (1.473 mag in the 



V band, see Raiteri et al. 20051 ). we can derive the cor- 
responding de-reddened flux density at 5450 A, F v = 
0.108 ± 0.017 mJy. Then X c L Xa = Aird 2 L vF v , where d L is 
the luminosity distance and A e the emitted wavelength. 
Assuming a flat cosmology with Hq = 71 kms -1 Mpc" 1 
and Qm = 0.27, d^ = 6141 Mpc. As a r esult, AL3000 ~ 
25.1 x 10 37 W. The IMcLure fc Dunlod (|2004 equation 



yields a best fit of i?BLR ~ 136 light days, with a range of 
possible sizes between 56 and 291 light days. 

4.1. Black hole mass 



Recently, iLiu et alJ (|2006l ) estimated the black hole mass 
for AO 0235+164 using several methods. In particular, 

4 From the WEBT archive. 



Mbh = G i?BLR ^blr> wnere Vblr is the Keplerian ve- 
locity of the line-emitting gas. The result wa s a black 



hole mass of 4.42 x 10° Mq. We notice that ILiu et al 



2006) derived the continuum flux from the Cohen et al 



19871 ) spectrum of October 1985, when the source was 
faint, but not at its minimum brightness level, and did 
not correct for Galactic and foreground galaxy absorption. 
But, chiefly, they implicitly assumed Vblr = FWHMMgn, 
which holds if the BLR gas has a disc-like geometry and 
if the angle between the l ine of sight and the disc axis is 
i = 30° (see discussion in IMcLure fc Dunlopl [20041 ) . This 
is an important point, since the black- hole mass depends 
on the squared gas velocity. But if the disc axis coincides 
with the jet axis, then the angle i should be very small for 
the str ongly beamed BL Lacs; for AO 0235+164 in par- 



ticular iFrev et al 



(2006) estimated i = 5° by analysing 
space VLBI maps during a faint state in 2001-2002. In this 
case, however, a disc-like geometry would imply that we 
should systematically see much narrower lines in blazars 
than in other AGNs (or that blazar black holes are system- 
atically much bigger). Hence, the assumption of randomly- 
oriented orbits for the BLR gas appears more appropriate, 
introducing a factor Vo/2 into the V-FWHM equation. 
For a BLR size of 136 light days and an FWHM Mg n ~ 
3300 kms -1 , as derived in the previous sections, we then 
obtain a black hole mass of ~ 2 x 10 8 Mq. On the other 
hand, with the less likely disc geometry and i ~ 5°, being 
V ~ FWHM/(2sini), we would obtain M BH ~ 1O 1O M , 
even if it can be intended as an upper limit. 

In conclusion, we stress that a black hole mass estimate 
based on this method suffers from many uncertainties, in 
particular those related to the unknown geometry of the 
BLR, so it should be considered with extreme caution. 



5. Microlensing 

To obtain significant amplification of the BLR by mi- 
crolensing, its size should be comparable to the projected 
radius of the Einstein ring of the lens on the source plane 
D s 9e, where 

4GM D ds 



c 2 D d D s 

is the angular radius of the Einstein ring, M the lens mass, 
and where -D s , D d , and D ds are the angular diameter dis- 
tances to the source, to the deflecto r and between the 
source and the deflector, respectively ( Einstein . 19361 ). 

With the cosmological assumptions mentioned in Sect. 
4, D s = 1632 Mpc, Ai = 1286 Mpc, and D ds = 621 Mpc, 
so that the Einstein radius of a 1 Af© star in the z = 
0.524 galaxy that microlenses emission from AO 0235+164 
is about 1.2 /xas. By comparison, a 136 light-day BLR 
subtends an angle of about 14 /ias at the redshift of AO 
0235+164. The resu l ting amplification is o nly about 1% 
(|Abaias et all l2002t ISchneider etail Il992j) and - 8% if 
the size is 56 light days. Moreover, the event would last 
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several years, since its duration is given by the time the 
lens takes to cross the projected radius of the Einstein ring 
on the lens plane: i var ~ where v is the transverse 

velocity of the star relative to the source-Ea rth line, whose 
typic al value is of the order of 10 3 kms _1 (jKavser et al 



1986). Significant amplification would be obtained in the 



case of a massive deflector, but this is expected to be a rare 
event and would imply even longer time scales. Actually, 
the most likely microdcflcctors in galaxies are thought to 
be massive compact halo objects (MACHOs) of subsolar 
mass. 

Hence, unless we have strongly overestimated the size 
of the BLR, we can conclude that we do not expect to 
detect variations in the Mgn broad emission line due to 
microlensing by a star. On the contrary, since the re- 
gion emitting the continuum radiation can be significantly 
smaller than the BLR, the possibility remains that mi- 
crolensing can affect the continuum flux. To have a typical 
2-3 month timescale event with a factor 2-3 flux ampli- 
fication, as observed during the big radio outbursts, we 
would need an ~ 10~ 4 Mq deflector acting on the radi- 
ation emitted from a very compact region of ~ 0.1 light 
dayfl 

More complicated scenarios arise when the optical 
depth of microlensing increases, so that one cannot con- 
sider isol ated stars, but rather their combined effect 



see 



e.g. iKavser et all Il986l ISchneider fe Weissl Il987 



Kavser et allll989tlTorres et alll2003h . The investigation 



of these scenarios goes beyond the aim of this paper, but 
we note that in the case of high optical depth, microlens- 
ing produces a sort of flickering rather than distinct, big 
events. Moreover, also de-amplification is possible under 
particular conditions. 

The fact that the foreground galaxy of AO 0235+164 
is not closely aligned with the blazar implies that lens 
candidates belong to the outer, less populated regions of 
the galaxy and dark halo, where the lensing optical depth 
is lower. 

6. Conclusions 

We have analysed five optical spectra of the BL Lac object 
AO 0235+164 taken at the VLT and TNG in 2003-2004. 
The source was in a faint state, and a broad Mg n emission 
line was clearly detected in all spectra. The flux of the 
line showed an overall variation by a factor 1.9, while the 
continuum flux density changed by a factor 4.3. 

The observed Mgn flux dependence on the contin- 
uum flux density may suggest that the line flux increases 
with source brightening but only until a certain satura- 
tion level, after which it may even decrease. A further 



Wc 



mention 



that, 



according 



to 

iGopal-Krishna fc Subramanianl (|l99ll h very fast variabil- 
ity events can be explained as due to microlensing if the 
radiation so urce is a superl u minal blob in the blazar jet. With 
this model, iRomero et al. I <| 19951 ) could explain the strong 
intraday radio variability of PKS 0537—441 observed in 1993 
with lens of 10~ 4 -10" 3 Mq. 



observing effort should be undertaken in order to check 
this trend. More likely, the line flux is independent of the 
optical continuum, a s found in the case of BL Lacertae 
( Corbett et all . [2000h . and it instead depends on another 
component. This can be identified with the one producing 
the UV-soft-X-ray bump in the source SED. 

In this case, the BLR-ionising continuum is most visi- 
ble in the optical spectral range when the source is in its 
faintest optical state. From the historical F-band mini- 
mum brightness, we thus estimated the size of the BLR 
by exploiting its dependen ce on the monochromatic lumi- 
nosity at 3000 A found bv lMcLure fe Dunlopl \2004) . The 
best-fit size is of the order of 4.5 light months, which has 
however to be taken as an upper limit. 

We investigated the possibility of explaining the flux 
variations we observed in the framework of a microlensing 
event due to a star in the foreground galaxy at z = 0.524. 
Although significant amplification of the continuum flux 
can also be obtained by subsolar-mass MACHO-like stars 
because of the very compact emission region, the size of 
the BLR we derived cannot account for a factor-2 line 
amplification unless we consider microlensing by a very 
massive star. Apart from the low probability of encoun- 
tering such an event, it would also imply variability time 
scales of the order of several years, in disagreement with 
the observations. 
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